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Top Mass Measurements from Jets and the Tevatron Top-Quark Mass 
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Theoretical issues are discussed for the measurement of the top-mass using jets, including perturbative and 
non-perturbative effects that relate experimental observables to the Lagrangian mass, and appropriate choices 
for mass schemes. Full account for these issues is given for e^e~ — > tt using a factorization theorem for event 
shapes for massive quarks. Implications for the Tevatron top-mass measurement are discussed. A mass-scheme, 
the "MSR-mass", is introduced which allows for a precise description of observables sensitive to scales 7? <C m, 
but at the same time does not introduce perturbative matching uncertainties in conversion to the MS mass. 



(N 
> 
(N 
(N 
(N 
p 

od 
O 
00 

O 



X 



1. Introduction 

The top-mass is a key parameter in the stan- 
dard model. For example, it plays an impor- 
tant role for analyzing electroweak precision con- 
straints, for predicting rare decays like B Xgj 
and K]^ — > TT^uiy, as well as for unraveling the 
Higgs sector in supersymmetric models. The lat- 
est measurement from the Tevatron [T], ml°^ — 
172.6 ± 0.8(stoi) ± l.l{syst) GeV, is of very high 
precision, and brings to mind several theoreti- 
cal questions. In what mass-scheme is the value 
quoted? Are measurements of the mass of a col- 
ored particle with a precision better than Aqcd 
possible using jets, or is this an irreducible un- 
certainty? How do perturbative and nonpertur- 
bative effects modify the relation between the 
experimental observable and the underlying La- 
grangian mass parameter? For the LHC a 
1 GeV systematic uncertainty was obtained from 
preliminary ATLAS studies [5] (where hurdles to 
go beyond this level include understanding the 
jet-energy scale to better than 1%). At this level, 
understanding the answers to the above questions 
is important. 

To achieve a high precision top-mass, measure- 
ments exploit kinematic information by consider- 
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ing leptons -I- jets, pp ttX {bqq){MD), or the 
associated dilepton or all-hadronic channels, and 
use Monte Carlo (MC) simulations to reconstruct 
observables sensitive to the top-quark four vec- 
tor, and hence to^ . One of the most sensitive ob- 
servables in this reconstruction are the invariant 
masses = (E.eaPf)' ^nd M| = (E.^.p^' 
of the jets and other decay products produced by 
the top and antitop respectively. Here the sets a 
and b depend on the jet-algorithm and cuts. In 
such analyses it might appear natural to think of 
the reconstructed top-mass as being associated to 
the pole- mass ?7t.|'°'°, since it is designed for sen- 



sitivity to the on-shell region pf 



However, 



as we will argue in section 5 below, it is not m^"^'' 
that is being measured by the Tevatron analyses. 

The relation of the pole-mass and any other La- 
grangian mass-scheme mt {R, can be expressed 
as a perturbative series, mf°'° ~ mt{R, fi) + 
5mt{R, fj,) where 



Smt{R,n) =R^Y1 



a-nk 



n=l k=0 



An 



with R a dimension- 1 scale intrinsic to the 
scheme, and ank finite numerical coefficients. 
Theoretically the pole-mass is a poor scheme 
choice since in QCD its definition is ambiguous 
by an amount of ©(Aqcd) (in perturbation the- 
ory this is referred to as the infrared renormalon 
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problem). Nice mass schemes avoid this prob- 
lem by a suitable choice for the a„fc's, and are 
known as short-distance mass schemes. Using 
observables expressed in terms of short-distance 
schemes the accuracy that the mass of a colored 
particle can be measured is not limited by Aqcd ■ 
This fact is important to obtain the ~ 40MeV 
uncertainties for current measurements of the b- 
quark mass [3J. For current top-measurements 
based on reconstruction there is another impor- 
tant restriction on viable short-distance schemes. 
In these kinematic-based analyses the top-quark 
decay is treated with a Breit-Wigner in the MC's, 
and only "top-resonance mass schemes" with a 
small i? ~ Ft are compatible with the Breit- 
Wigner form [4]. This restriction rules out di- 
rectly measuring the MS mass, which has a much 
larger R = mt. 

To explore these issues in detail it is useful to 
consider a case with full analytic control, namely 
jets produced in e+e~ — > tt with cm. energy 
Q ^ nit. Here the jet-invariant masses Mj^ 
sum over particles in the top and antitop hemi- 
sphere defined with respect to the thrust axis, and 
are examples of event shape variables for massive 
quarks. The resonance region jAfj j — mt| <C rrit is 
most sensitive to the top-mass. The appropriate 
factorization theorem was derived in Ref. [4], 

^^^1^ =aoi/Q(g,M™)i^™(^,|,M™,MA) 

X / d£+drB+{St ~9^,Sm, Ft, ^^a, Mr) (2) 
J m 

xB_{si~^,6m, Tt, MA, fir)S{£+,r,SA, A, ma), 
m 

where St = {Mf — m?)/m and si = {Mf — 
m^)/m. Eq. ([2]) is valid to all orders in as and 
has power corrections of 0{mas{m)/Q,m^ /Q^, 
Fj/m, St^f/m). Here Hq and Hm are pertur- 
bative coefficients describing the hard scales Q 
and m, as well as the summation of large log- 
arithms from Q ^ TO ^ Mr, where Mr ^ 
Ft + Q Aqcd /to -I- Sf. The B± are perturbatively 
calculable heavy-quark jet functions, which de- 
scribe the evolution and subsequent decay of the 
top/antitop quark to jets, and include a resum- 
mation of large logs between Mr and the soft scale 



Ma ^ Aqcd + rnTt/Q + ms/Q. A complete re- 
summation of large logs at next-to-leading or- 
der (NLL) was carried out in Ref. [5], where the 
reader can find further details. In Eq. ^ m 
and Sm are specified in a top-resonance mass- 
scheme, of which an example known as the jet- 
mass is discussed below in section 2. Finally, S is 
the hemisphere soft-function describing the soft- 
radiation between jets. S is perturbatively cal- 
culable for £^ ^ Aqcd and non-perturbative for 
£^ ~ Aqcd- The parameter A is a gap in the soft 
function [6], while 6 A indicates that the scheme 
for this gap is defined so that both perturbative 
and nonperturbative contributions are included 
without inducing a leading renormalon ambigu- 
ity, as discussed below in section 3. 

The result in Eq. ([2|) relates an experimentally 
measurable hadronic observable, (Pa/dM^dM^ , 
to the short distance mass m. As an example, 
consider the peak position of the invariant mass 
distribution Mt^i = MP"^"". Schematically this 
relation has the form 

MP-'^ = m + Ft(a,+a2 + ...) + ^^^, (3) 

TO 

where the perturbative shifts ^ Tt{as + + . . .) 
can be calculated from the jet-functions, B±, and 
the nonperturbative shift ~ QAqcd/to is com- 
pletely determined by the soft-function S. Hav- 
ing a systematic separation of the perturbative 
and non-perturbative shifts into B± and S is im- 
portant for a precision determination of the La- 
grangian mass to. In particular the value of the 
terms in the perturbative series in Eq. ([3]) are re- 
lated to the choice of mass-scheme for to, and are 
modified by different scheme choices. Thus, con- 
trol over the mass-scheme requires control over 
these corrections. Note that the non-perturbative 
shift is always positive for hemisphere masses, 
and quite sizeable, ~ 1-2 GeV for typical values 
of Q. In sections 2 and 3 below we discuss the 
functions responsible for these perturbative and 
nonperturbative shifts in more detail. Then in 
section 4 we quantitatively analyze the relation 
between MP"""^ and TO using Eq. Finally, in 
section 5 we return to implications for the Teva- 
tron mass, and discuss a very useful mass-scheme 
for measuring both top and bottom heavy-quark 
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masses, the MSR-scheme. 

2. Heavy-Quark Jet Function 

Since s <^ m the jet function can be formulated 
in the heavy-quark hmit with HQET, 



C 



HQET 



h„ i iv ■ D 



5m + ^ ) hr. 



(4) 



Here Sm = mP°'° — m encodes the mass-scheme 
choice, and Ft is the total width of the top-quark 
which suffices since our observable is inclusive in 
the decay products. The heavy-quark jet function 
is defined as the imaginary part of a hy prop- 
agator connected by a light-like Wilson line [4], 
B{s, 6m, Tt, /i) = Im[B(s, Sm, Fj, /i)] where 



B{s,Sm,rt, fi) = 



(5) 



{0\TKiO)Wn{0)W^{x)K{x)\0). 



At tree level in as we have a Breit-Wigner, B = 
Ft/[(7rm)(s2-HF^)]. From Eq. (g]) the jet function 
obeys the shift identity B{s, Sm,Tf , fj.) ~ B{s — 
25m-\-iTt, 0, 0, y) = B{s—2Sm+iTt, n). It also has 
an RGE /Ltd/(i/ii3(s, /i) = J ds'jB{s~s' , fi)B{s' , ^) 
where the anomalous dimension has a plus- 
function term with the universal cusp anoma- 
lous dimension and a (5-function term, 7s (s, /i) — 
-2F™''P[a,]//i[/i0(J)/s]+ +7B[a,]5(s). Summing 
logarithms from /ir to /iA gives 

B{s,Sm,Tt,^J,A,^J'^) 

Ids' UB{s-s',fiA, fJ-r)B{s', Sm,Tt, fir) , (6) 



which is the function B — B 



+ 



B. 



appearmg 



in Eq. Eq. ([5]) sums all logs that can affect 
the shape of d^a/dM^dMf [S]. There is some 
freedom in combining fixed order results for B at 
Hr with the RGE in Ub- Wc define a counting 
that elevates the importance of fixed order results 

LL: 1-loop Fcusp, tree-level matching ; (7) 
NLL: 2-loop Fcusp, 1-loop 7s and matching; 
NNLL: 3-loop Fcusp, 2-loop 7s and matching. 

One-loop results for B and 7s were computed in 
Ref. [S], and two-loop results in Ref. [7]. The 
three-loop result for Fcusp is known from Ref. [5] , 
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Figure 1. Heavy quark jet-function up to NNLL 
order (top). Peak positions in the pole and jet- 
mass schemes (bottom). Plots from [7]. 

so B is known at NNLL. Results for B in mo- 
ment space have also been computed [1019] . A 
different heavy quark jet function also occurs in 
the massive quark form factor, whose anomalous 
dimension shares several coefficients in common 
with 7s [5]. 

In order to implement a short-distance mass- 
scheme and subtract the infrared contributions 
related to the pole- mass renormalon we must sys- 
tematically expand in dm to the same level in as 
that we determine B itself. Writing 



Sm 



5 (a*) 



Smi 



Sm 



2}, (8) 



with R Tt, the remaining freedom in specifying 
the mass-scheme at two-loop order corresponds to 
defining Smi^2- The jet-function itself can be used 
to define a jet-mass scheme, m = mj(R,fj,) [1]. 
It is convenient to use the position space jet- 
function B(y,fi) so that mj has a consistent /z- 
anomalous dimension 17 , and we define 



Smj 



_ Re''' d\nB{y,n) 
^ ~ d\n[iy) 



iyefE=l/R 



(9) 
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This definition for the jet mass is possible for any 
value of R, and the anomalous dimension in R 
can also be consistently derived [TT] , 

In Fig. [T] (top panel) we show results for Eq. ^ 
using R — 0.8 GeV, the jet-mass with reference 
value mj(i?, 2.0 GeV) = 172 GeV, ha = 1 GeV, 
and three curves at each of LL, NLL, and NNLL 
order corresponding to /zr = 3.3, 5.0, 7.5 GeV. 
Convergence by NNLL is evident. In the bot- 
tom panel of Fig. [T] we compare the pole and 
jet-mass schemes and show that very good con- 
vergence for the peak position is achieved in the 
jet-mass scheme. 
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Figure 2. Diagonal hemisphere soft function up 
to NNLO order in the 5 A = scheme and scheme 
from Eq. (HH). Plot from [13]. 



3. Hemisphere Soft Function 

The hemisphere soft function is defined by a 
matrix element of Wilson lines 

5hemi(^+, r,/i) - -^5]^(^+-fcr)'5(r-fc7'') 



x(o|rft(o)r„(o)|x,)(x,|y,t(o)rl(o)|o), (lo) 

where Y,l{0) = Pexp{ig J^dsn-A{ns)) and 
is similar but in the 3 representation. Here fc+'^ 
is the sum of all plus-momenta of particles in Xs 
in hemisphere a, and k'^^ is the sum of all minus- 
momenta for particles in hemisphere b. Thus the 
(5-functions in Eq. (jlOp are a refiection of the pre- 
scription for assigning soft-radiation to the invari- 
ant masses Mt and M^. This soft function is uni- 
versal [4], the same function appears in the fac- 
torization theorem for the hemisphere invariant 
mass distribution as well as the thrust and heavy 
jet mass distributions of massless jets [12j . 

S{£'^,£^,fi) has both perturbative and non- 
perturbative components, and a convenient way 
to account for this is [6J 

sie+,e-,SA,A,fiA) = Jm+'m-' (u) 

xSP''"' {£+-£+', r-r\ 5 A, ha)S'"°'^{£+\ £-', A) . 

Here 5"^^'* refers to a partonic computation of 
Eq. pU]) and encodes the proper ^a dependence 
and large £^ behavior into S. 5'™°'^ is a model 
for the non-perturbative £^ ^ Aqcd region of 
S. A gap 9{£+' - A)9{£-' - A) is contained 



in 5*™°^ where A is defined in a scheme spec- 
ified by (5A(i?',/i). The partonic soft-function 
obeys S'p"*(^+, JA, ^) = SP'''^{£+ - SAJ- - 
<5A,0,^) = S'P"*(£+ - SAJ- - 6A,fi). To avoid 
a O(Aqcd) renormalon (which has nothing to do 
with the ©(Aqcd) renormalon of quark masses), 
we must systematically expand in SA to the same 
level in as as we expand S^'^'^^ [B] . We write 



SA^R'e^^\^^6Ai 

I TT 



6A^ 



(12) 



Removing the renormalon with SA avoids having 
to make large changes to 5'™°'^ each time an ad- 
ditional order in as is included in S^^'^^. Much 
like Eq. ^ we can define ^A with the position 
space soft function [T31, which yields a scheme 
for A{R, fi) with a consistent RGE in fi and R 



6 A = R'e'^" 



d\n{ixi) 



ie"'Exi^2=l/R' 



(13) 



Results for S'^^^^{£'^ ,£~ ,SA, ^) are known at one- 
loop order [5] and two-loop order [T3] , and for the 
thrust soft-function in Refs [14]. 

In Fig. [2] we compare results for Eq. ([TT]) at 
LO (black), NLO (blue), and NNLO (red), in the 
SA = scheme (dotted & dot-dashed) and in 
the A-scheme (dashed & solid). For S''^°^ we use 
the 2d-exponential model of Ref. [15] with the 
addition of the gap A. The negative dips at small 
£ that are present in the SA — scheme, are 
removed in the renormalon free A-scheme, and a 
stable result is obtained for the peak-position of 
the soft function. 
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Figure 3. Diagonal invariant mass distribution at 
LL and NLL order (top). Peak position of the 
NLL distribution for different soft function mod- 
els versus Q/m (bottom). Vertical black lines in- 
dicate theory uncertainties and the colored bands 
linear extrapolations to Q/m ^ 0. Plots from [5]. 



4. NLL e+c 



tt Cross-Section 



We now proceed to put together all the ingre- 
dients in Eq. ([2]) at NLL order [5J (since the full 
NNLL results have not yet been published). The 
cross-section depends on /xg ~ Q, [im ^ rn, /zr, 
and fiA, and the dependence on these four renor- 
malization scales is reduced from LL to NLL or- 
der. The largest variation occurs for /ir and /iA, 
but is highly correlated. It becomes significantly 
smaller if we vary fixing ^r/z^A = Q/rn, which is 
consistent with the scaling given below Eq. 
Figl3] (top panel) shows the invariant mass dis- 
tribution for Mt = Ml where the three curves at 
each order have /zr = 3.3, 5.0, 7.5 GeV and the 
same reference top- mass of 172 GeV was used as 
in Fig. [2 The O(Aqcd) renormalons associated 
to the pole-mass and the partonic soft-function, 



both would cause bigger shifts to the location 
of the invariant mass peak than the difference 
(~ 0.15 GeV) between LL and NLL order shown 
in Fig. [3) Thus using mass and gap schemes 
that avoid these renormalons was important to 
obtain perturbative stability. The largest correc- 
tion to ^'^ — m is due to the soft radiation, 
which causes the majority of the difference be- 
tween the dashed and solid vertical lines of FigO 
(top panel). The advertised shift, oc QAqcd/w, 
is demonstrated in Fig. [3] (bottom) where we show 
that linear dependence on Q occurs for six differ- 
ent models for S"™"*^. For all models the extrap- 
olations to Q/m converge to the underlying 
Lagrangian mass value m. In an e^e~ environ- 
ment the correction Mf^^^ — m can be determined 
a) from measurements of S for massless jets that 
together with perturbative computations deter- 
mine the shift as in Fig. [3] (top), or b) from mea- 
surements at different Q's with a linear extrapo- 
lation that removes the soft-radiation effect as in 
Fig. [3] (bottom). 

5. Implications for the Top Mass in MC 
programs and the Tevatron Top Mass 

The situation at high energy hadron colliders 
is somewhat different from e^e~ collisions due 
to the more complicated initial state and because 
the jet algorithms needed for the reconstruction of 
the invariant mass distribution are implemented 
in terms of MC programs rather than a factoriza- 
tion theorem. Nevertheless, the principles upon 
which final state interactions in MC's for tops 
are based, have close analogues with elements in 
the factorization theorem These analogies 

make it obvious that the top quark masses con- 
tained in MC programs, which are measured at 
the Tevatron and LHC, are very similar to the jet 
mass definition we discussed above. Moreover, 
the Tevatron top mass measurements are statis- 
tically dominated by the reconstructed invariant 
mass distribution in the peak region, which al- 
lows more reliable statements to be made about 
the scheme of the Tevatron top mass. 

The final state shower in MC's describes the 
perturbative aspects of coUinear and soft radia- 
tion of partons. Starting at transverse scales of 
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order the momentum transfer of the primary par- 
tons the parton shower evolves the system down 
to the shower cutoff scale i?sc which is typically 
in the 1 GeV range. For the top quark the width 
Tt ~ 1.5 GeV provides an additional natural 
shower cutoff for radiation off the top quark. In 
the factorization theorem ([2]) the same physics is 
described by the renormalization group running 
of the factors Hq and iJ„j , the partonic contribu- 
tions in the soft function and, in particular, the 
jet functions B±. As far as the question of the 
mass scheme implemented in MC's is concerned, 
the analogue of the shower cutoff Rsc is the modi- 
fication of the jet functions B± due to the residual 
mass term Sm, which subtracts low-energy fluc- 
tuations of the jet functions, absorbing them into 
the mass definition. It is therefore the implemen- 
tation of the parton shower for the top quark and 
the size of the shower cutoff Rsc which determines 
the top mass definition used in the MC's. 

Another important ingredient of MC's is the 
description of nonperturbative effects in the 
hadronization process through models that de- 
pend on many parameters. These are fixed from 
reference processes. For the factorization theo- 
rem ^ the analogue of the hadronization models 
is the soft function which can be determined from 
event-shapes involving light quark jets. We note 
that at hadron colliders the treatment of nonper- 
turbative effects is more involved due to the par- 
tonic initial state, and does not necessarily lead 
to a positive shift of the peak position. 

For hadron collisions at the Tevatron and the 
LHC there are additional complications due to for 
example underlying events, more involved com- 
binatorial background and initial state showers. 
However, if such effects are modeled correctly in 
the MC's they do not affect the correspondences 
mentioned above. We therefore conclude that the 
top quark mass that is implemented in MC's has 
the property 



Rs 



(14) 



with Rsc ~ Ff '--^ 1 GeV. The exact coefficients 
in the series on the RHS depend on how the 
parton shower is implemented and are currently 
unknown. However, it is reasonable to assume 




150 R 



Figure 4. Converting the Tevatron top mass into 
the MS scheme using the MSR scheme. 

that they are of order one. Numerically, the MC 
masses are therefore quite close to the jet mass 
scheme defined in ([S]) for which R ^ Rsc or to 
any other short-distance mass scheme that falls 
in the category of ^ with R ~ i?sc- 

A good way to illustrate the numerical size of 
the scheme uncertainty due to the ignorance of 
the coefficients in (fT4|) is the associated uncer- 
tainty in the MS top mass Wit(jnt)- This uncer- 
tainty is essential for fits of electroweak precision 
observables and the resulting indirect Higgs mass 
bounds. To avoid large logarithms in the pertur- 
bative relation of the MC or the jet mass (with 
i? ~ 1 GeV) and the MS mass (with i? ~ mt), one 
can use the evolution equation in R from Ref. [llj . 
Conversions involving the MS mass rnt{rnt) are 
particularly simple when the R-evolution is car- 
ried out with the MSR- mass scheme m^^^(i?), 
which is a short-distance mass scheme that also 
falls in the category of ^ and has in addition 
the property ^^(mf ^^) = rnt{rrit). The MSR 
mass definition is obtained from the MS-pole se- 
ries by a simple replacement rule [TT]. We can 
identify the MSR mass at a low scale R in the 
1 GeV range with the mass that is contained in 
MC programs. 



mf^iRsc) = ™r^(3i^ GeV) 



(15) 



The variation of R between 1 and 9 GeV param- 
eterizes our remaining ignorance. 

We now identify mf^^Rsc) with the Teva- 
tron measurement m}^'^ = 172.6 ± 0.8(siai) ± 
l.l(s?;si) GeV 1 . The result for rnt{mt) is il- 



7 



lustrated in Fig. U) The vertical error bars on 
the Tevatron mass are the experimental uncer- 
tainties and the horizontal error bars reflect the 
current scheme uncertainty due to Eq. (|15p . The 
three red lines show the 3-loop R-evolution of the 
scheme uncertainty with the MSR mass. The in- 
tersection of the red lines with the black line that 
shows the equation m^^^ {R) = R, gives the MS 
top mass mt(mt) = 163.0±1.3+^:^ GeV [H]. Here 
the first error is the combined experimental one, 
and the second error is from the scheme uncer- 
tainty in the Tevatron mass. 

The above correspondence between the shower 
MC and the factorization theorem is compelling 
for large pt events where top quarks are fast. A 
potential concern is that tops at the Tevatron 
are produced with predominantly small px and 
are therefore slow. However given that the MC 
framework applies to this situation, the concept 
of the top-mass intrinsic to the MC should be 
independent of whether the MC is applied to en- 
ergetic tops or to soft tops, and hence the corre- 
spondence we have discussed should be applicable 
for the Tevatron as well. Further study of this is 
warranted. 
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